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Abstract 

Magnetite (Fe304), an archetypal transition metal oxide, has been used for thousands of years. 



from lodestones in primitive compassesyj] to a candidate material for magnetoelectronic devices.|l2t] 
In 1939 Verweyl^] found that bulk magnetite undergoes a transition at Ty ~ 120 K from a high 
temperature "bad metal" conducting phase to a low-temperature insulating phase. He suggestedyfl 



that high temperature conduction is via the fluctuating and correlated valences of the octahedral iron 
atoms, and that the transition is the onset of charge ordering upon cooling. The Verwey transition 
mechanism and the question of charge ordering remain highly controversial.] 5, ^ 7, 8, ^ 10, 11] 
Here we show that magnetite nanocrystals and single-crystal thin films exhibit an electrically driven 
phase transition below the Verwey temperature. The signature of this transition is the onset of sharp 
conductance switching in high electric fields, hysteretic in voltage. We demonstrate that this transition 
is not due to local heating, but instead is due to the breakdown of the correlated insulating state when 
driven out of equilibrium by electrical bias. We anticipate that further studies of this newly observed 
transition and its low-temperature conducting phase will shed light on how charge ordering and 
vibrational degrees of freedom determine the ground state of this important compound. 
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Strongly correlated electronic materials can exhibit dramatic electronic properties {e.g., high 
temperature superconductivity, metal-insulator transitions, and charge ordering) not present in 
simple systems with weaker electron-electron interactions. Such rich electronic phenomenology 
can result when electron-electron interactions, electron-phonon interactions, and electronic band- 
width are all of similar magnitude, as in magnetite. 11 1211 VerweyyD found nearly seven decades ago 
that bulk magnetite, while moderately conductive at room temperature, undergoes a transition to 
a more insulating state below what is now called the Verwey temperature, Ty ~ 120 K. Similar 
transitions are known in a number of materials. ifisi 



1411 Above Ty, Fe304 has an inverse- spinel 
structure of the form AB2O4, with tetrahedrally coordinated A sites occupied by Fe'^+ and oc- 
trahedrally coordinated B sites of mixed valence, equally occupied by irons with formal +3 and 
+2 charges. Conduction at high temperatures has long been thought to be through fluctuating 
valences of the B sites, with the transition corresponding to some kind of B site charge ordering 
as T decreases; concurrent is a first-order structural phase transition to an orthorhombic unit cell. 
This explanation remains controversial,! 5, ^ with experiments showing some charge dispropor- 
tion or charge order (C0),||7L 8|land others implying that the structural degrees of freedom drive 
the change in conductivity. [|9l[lOD Recent theoretical progress has been made in understanding the 
complex interplay of charge and structural degrees of freedom! 15, 1^, including a complete pic- 
ture of the transition mechanism 111 ill with strongly correlated 3d Fe electrons acting to amplify 
electron-phonon couplings. Testing these ideas experimentally is of much interest. 

In this Letter we report electronic transport measurements in magnetite at the nanoscale on 
both nanocrystals and single-crystal epitaxial thin films. Both types of devices exhibit striking 
electrically-driven hysteretic switching of the electronic conductance once sample temperatures 
are reduced below Ty. The data clearly show that the transition is not the result of local heat- 
ing above Ty, but instead is an electrically-driven breakdown of the insulating state. We dis- 
cuss possible explanations for this switching in the context of the general Verwey transition prob- 
lem. While qualitatively similar resistive switching has been observed in other correlated oxide 
systems n 171 1 1 811 . the phenomenon in Fe304 is a bulk effect with a mechanism distinct from these. 

Two-terminal devices for applying voltages and measuring conduction at the nanoscale have 
been fabricated (see Methods) incorporating both Fe304 nanocrystalsjl^ (10-20 nm in diam- 
eter with oleic acid coating) and single-crystal thin films (40-60 nm thick). 1I20I1 Devices were 
measured in both a variable temperature vacuum probe station and a ^He cryostat with magnet. 
Current- voltage characteristics have been measured with both a semiconductor parameter analyzer 
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FIG. 1: Hysteretic conductance switching below Ty. (a) Current- voltage characteristics at various temper- 
atures for a device based on 10 nm magnetite nanocrystals. Arrows indicate the direction of the hysteresis. 
Inset: Zero-bias resistance, R{T). (b) Analogous data for a device based on a 50 nm-thick MBE-grown 
magnetite film. The nominal interelectrode gap was planned to be 100 nm, but at its narrowest was ap- 
proximately 10 nm. Inset: Two-terminal resistance as function of channel length for another set of devices 
fabricated on another piece of the same film. 



and directly using voltage sources and current amplifiers, with differential conductance computed 
numerically. 

Figure la shows I — V characteristics of a nanocrystal device at selected temperatures. When 
cooling, zero-bias conductance decreases monotonically until T ^ Ty. Below Ty, the I — V 
characteristics show sharp switching between a low bias insulating state and a high bias state with 
much higher differential conductance dI/dV{V) (close to dI/dV{V = 0,T = 300 K)), with 
dramatic hysteresis as a function of voltage sweep direction. The switching threshold voltages 
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FIG. 2: Differential conductance plots of the switching, (a) dl/dV vs. V for the nanocrystal device from 
Fig. la at the same temperatures, (b) dl/dV vs. V for the thin film device from Fig. lb. Both plots have 
logarithmic dl/dV axes to better show the lowest temperature data. 



increase in magnitude as T is decreased. 

Dozens of nanocrystal devices were measured and only those with 300 K resistances below 
10 kfi showed the switching, with higher resistance devices having higher switching threshold 
voltages. Resistances decrease by some three orders of magnitude with vacuum annealing at 



673 K, likely because of oleic acid decomposition. [|2 IL 122(1 The temperature dependence (Fig. la, 



inset) of the zero-bias resistance, R{T), has no step at Ty, showing that R{T) remains dominated 
by contact effects. 

Qualitatively identical conduction is apparent in the thin film devices, as shown in Fig. lb. 
Contact resistances are also important in these structures, as demonstrated by examining R at low 
bias (< 100 mV) as a function of channel length, L, as shown in the inset for one set of devices. 
Extrapolating back to L ^ 0, the contact resistance, -Rc, at 300 K is 390 fi, while the 50 nm 
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thick channel of width 20 /um contributes 27.2 fi/micron, implying (based on channel geometry) a 
magnetite resistivity of 2.9 mfi-cm, somewhat below bulk expectations. Further investigations are 
seeking to understand and minimize Re- Analysis of R{L) at lower temperatures shows that Rc 
increases significantly as T is decreased, exceeding 80 kTl by 80 K. This complicates the analysis 
of the switching, since some of the total V is dropped across Rc rather than directly within the 
Fe304; further, the contacts may not be Ohmic near the switching threshold. We return to this 
issue later. 

The transitions in all devices are extremely sharp, with widths less than 50 /iV, though in 
repeated sweeps at a fixed temperature, there is sweep-to-sweep variability of a few mV in switch- 
ing thresholds. Using the substrate as a gate electrode, no discernable gate modulation was seen 
in nanocrystal devices for gate biases between -80 V and -1-80 V; this suggests that nanocrystal 
charging effects do not dominate. Switching characteristics were independent of magnetic field 
perpendicular to the sample surface up to 9 T, showing no large coupling between magnetization 
and the transition. 

Differential conductance traces (Figure 2) show the transition even more dramatically. In the 
high conductance state, dl/dV is relatively temperature independent. As T is decreased, a clear 
zero-bias suppression develops, deepening into a hard gap when T < Ty. In the nanocrystal 
data there are indications (in (Pi /dV'^{V)) of gap formation even at 150 K. We note that Ty in 
nanocrystals could be elevated, since nanocrystals have large surface-to-volume ratios and the 
transition temperature of the magnetite surface is known to be higher than in the bulk. [23] 

Several lines of evidence indicate that these sharp conductance transitions are not the result of 



loca 



Joule heating (as in macroscopic samples of Fe304 [|24L 



2511 and in the Mott insulator VO2 



2711 '). but rather are electrically driven. In the worst-case scenario, all of the I x V Joule heat- 



ing power is dissipated within the magnetite, and inhomogeneous dissipation (e.g., filamentary 
conduction through a locally heated path) can complicate the analysis. The local steady-state tem- 
perature is determined by the power dissipated and the thermal path. Thermally driven switching 
would then correspond to raising the local temperature above Ty. At a fixed cryostat temperature 
an improved thermal path would imply that more power dissipation would be required for a given 
local temperature rise. Similarly, for a fixed thermal path, the necessary dissipated power for ther- 
mal switching would approach zero as T — > Ty. Furthermore, at a given cryostat temperature 
thermally-driven switching would imply that the power dissipated at the low-to-high conductance 
transition (needed to raise the local temperature to Ty) should be close to that at the high-to-low 



5 




FIG. 3: Power required to switch from the insulating into the more conducting state as a function of tem- 
perature, for a device based on ~20 nm diameter nanocrystals. Inset: Hysteresis loop in the conduction of 
the same device at 80 K, showing essentially no change in switching characteristics as the voltage sweep 
rate is varied over two orders of magnitude. 



conductance transition. 

The theraial conductivity, k, of magnetite is dominated by phonons in this temperature range, 
and limited by phonon-electron scattering, pSQ even when T > Ty. As a result, k increases as T is 
decreased through and below Ty, and the material's thermal coupling to the cryostat improves as 
T is reduced. In all devices showing switching, the electrical power required to switch from low to 
high conductance decreases with decreasing T, with Fig. 3 showing one example. This is precisely 
the opposite of what one would expect from thermally-driven switching. Similarly, in all devices 
the power dissipated at switching does not approach zero as T — > Ty, again inconsistent with 
thermally-driven switching. Furthermore, at a given T the power dissipated just before V is swept 
back down through the high-to-low conductance threshold significantly exceeds that dissipated 
at the low-to-high point in many devices, including those in Fig. 1, inconsistent with thermal 
switching expectations. Finally, nanocrystal and thin film devices show quantitatively similar 
switching properties and trends with temperature, despite what would be expected to be very 
different thermal paths. These switching characteristics are also qualitatively very different from 
those in known inhomogeneous Joule heating. [|24l|25|]. These facts rule out local heating through 
the Verwey transition as the cause of the conductance switching. 

Figure 3 (inset) shows details of hysteresis loops on a nanocrystal device comparing different 
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Channel Length L [nm] 

FIG. 4: Switching voltages in a series of film devices as a function of channel length at several temperatures. 
The hnear variation with L strongly implies that for each temperature there is a characteristic electric field 
required for switching. The non-zero intercepts of the trend lines indicate that some device-dependent 
threshold voltage must be exceeded for switching even when L = 0, suggestive of contact effects. 



voltage sweep rates. The loop shape and switching voltages are unchanged to within the precision 
of the data collection as voltage sweep rates are varied from around 0.7 V/s up to 70 V/s. This 
indicates that the switching process is relatively rapid. Further studies will examine the intrinsic 
switching speed. 

The observed conductance transition appears to be driven electrically. Figure 4 is a plot of 
the low-to-high conductance switching voltage as a function of L in a series of film devices for 
several temperatures. The linear dependence implies that the transition is driven by electric field 
itself, rather than by the absolute magnitude of the voltage or the current density. The fact that the 
voltage extrapolates to a nonzero value at L = is likely a contact resistance effect. Minimizing 
and better understanding the contact resistance will allow the determination of the electric field 
distribution within the channel. 

The length scaling of the transition voltage also demonstrates that this is a bulk effect. The 
contacts in all of these devices are identical, so any change in switching properties must result from 
the magnetite channel. This is in contrast to the resistive switching in Pro.yCao.sMnOs (PCMO) 
that is ascribed to a change in contact resistance due to occupation of interfacial states! isll. 

The field-driven conductance transition may give insights into the equilibrium Verwey tran- 



sition. This switching may be useful in testing recent calculations 
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about the role of 
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strongly correlated B-site Fe 3d electrons and their coupling to phonons in the Verwey transition 
mechanism. It is interesting to ask, to what degree is the field-driven electronic transition coupled 
to the local structure? It is greatly desirable to perform local probes of the magnetite structure 
(via x-ray or electron-diffraction techniques or scanned probe microscopy) in situ in the channel 
of biased devices, to see if the coherence between structural symmetry changes and the formation 
of a gap in the electronic spectrum is broken under these noequilibrium conditions. This is a sig- 
nificant experimental challenge. Similarly, local Raman spectroscopy of devices under bias could 
reveal field-induced changes in phonon modes and electron-phonon couplings, and single-crystal 
thin films permit the application of bias along well-defined crystallographic directions relevant to 
structural symmetry changes at Ty. We do note that qualitatively identical switching occurs in 
nanocrystal devices as in strained thin films strongly coupled to rigid MgO substrates. This sug- 
gests that elastic constraints on scales much larger than the unit cell have relatively little influence 
on the observed switching. 

It is also possible that the nonequilibrium carrier distribution contributes to destabilizing the 
insulating state. In the presence of a strong electric field a carrier can gain significant energy even 
in a single hopping step, even though carrier relaxation times are very short. A rough estimate of 
the average critical £'-field for switching at 80 K is 10'' V/m, from the slope of the line in Fig. 4. 
The high temperature cubic unit cell is 0.84 nm on a side, meaning that a carrier traversing one 
cell would gain approximately 8.4 meV, comparable to k^Ty ^ 10.3 meV. Conductance switching 
at such high fields may require consideration of such nonequilibrium carrier dynamics. 

The presence of multiple switching transitions in individual nanocrystal and film devices also 
bears further study. The suggested charge order may melt inhomogeneously, with portions of 
the channel having different switching thresholds. There could also be charge-ordered interme- 
diate states between the insulating regime and the most conducting regime. [5] Again, optical 
measurements n 1 211 with sufficient spatial resolution could address these possibilities. Through 
improved metal/magnetite contacts and further study, it should be possible to unravel the pre- 
cise nature of this nonequilibrium transition, and its relationship to the equilibrium, bulk Verwey 
transition. 
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I. METHODS 



Magnetite nanocrystals were prepared via solution-phase decomposition of iron carboxylate 
salts . [1 1911 The nanocrystals have been characterized by transmission electron microscopy (TEM), 
x-ray diffraction, and infrared and Raman spectroscopy, as discussed in Supplemental Material. 
As synthesized the nanocrystals are protected by weakly bound oleic acid ligands; these ligands 
allow the suspension of the nanocrystals in organic solvents, but act as electrically insulating layers 
that must be largely removed for effective electronic transport measurements. 

Two-step electron beam lithography and e-beam evaporation (1 nm Ti, 15 nm Au) were used 
to pattern closely spaced source and drain electrode pairs onto degenerately n-doped silicon sub- 
strates coated with 200 nm of thermally grown Si02. Interelectrode gaps (channel lengths) ranged 
from zero to tens of nm, with a 10 fim wide channel region. Nanocrystals were spin-coated from 
hexane solutions to form slightly more than one densely packed monolayer of nanocrystals over 
the channel region. Samples were then baked at 673 K in vacuum for 1 hr to remove as much of 
the oleic acid as possible. In one set of samples, a second layer of particles was added followed 
by a second round of baking. 

The other class of devices are based on epitaxial magnetite films 50 nm thick grown by oxygen- 
plasma-assisted molecular beam epitaxy (MBE) on (100) MgO single-crystal substrates. Details 
of the growth process have been reported elsewhere. [20] Single-step e-beam lithography and e- 
beam evaporation were used to pattern Au (no Ti adhesion layer) source and drain electrodes 
defining a channel length ranging from tens of nm to hundreds of nm, and a channel width of 
20 fim. The interelectrode conduction is dominated by the channel region due to this geometry. 
No annealing was performed following electrode deposition. 



[1] Mills, A. A. The lodestone: history, physics, and formation. Ann. Science 61, 273-319 (2004). 

[2] Coey, J. M. D. & Chien, C. L. Half-metallic ferromagnetic oxides. MRS Bulletin 28, 720-724 (2003). 

[3] Verwey, E. J. W. Electronic conduction in magnetite (Fe304) and its transition point at low tempera- 
tures. Nature 144, 327-328 (1939). 

[4] Verwey, E. J. W. & Haayman, P. W. Electronic conductivity and transition point of magnetite. Physica 
8, 979-987 (1941). 

[5] Walz, F. The Verwey transition - a topical review. /. Phys.: Condens. Matter 14, R285-R340 (2002). 

9 



[6] Garcia, J. & SuMas, G. The Verwey transition - a new perspective. J. Phys.:Condens. Matter 16, 
R145-R178 (2004). 

[7] Huang, D. L. et al. Charge-orbital ordering and Verwey transition in magnetite measured by resonant 

soft x-ray scattering. Phys. Rev. Lett. 96, 096401 (2006). 
[8] Nazarenko, E. et al. Resonant x-ray diffraction studies on the charge ordering in magnetite. Phys. Rev. 

Lett. 97, 056403 (2007). 

[9] Subfas, G. et al. Magnetite, a model system for mixed- valence oxides, does not show charge ordering. 
Phys. Rev Lett. 93, 156408 (2004). 
[10] Rozenberg, G. K. et al. Origin of the Verwey transition in magnetite. Phys. Rev. Lett. 96, 045705 
(2006). 

[1 1] Piekarz, R, ParUnksi, K. & Oles, A. M. Mechanism of the Verwey transition in magnetite. Phys. Rev. 
Lett. 97, 156402 (2006). 

[12] Gasparov, L. V. et al. Infrared and Raman studies of the Verwey transition in magnetite. Phys. Rev. B 
62, 7939-7944 (2000). 

[13] Imada, M., Fujimori, A. & Tokura, Y. Metal-insulator transitions. Rev. Mod. Phys. 70, 1039-1263 
(1998). 

[14] Coey, M. Condensed-matter physics: Charge-ordering in oxides. Nature 430, 155-157 (2004). 

[15] Leonov, 1., Yaresko, A. N., Antonov, V. N. & Anisimov, V. 1. Electronic structure of charge-ordered 

Fe304 from calculated optical, magneto-optical Kerr effect, and O K-edge x-ray absorption spectra. 

Phys. Rev. B 74, 165117 (2006). 
[16] Pinto, H. P. & ElUot, S. D. Mechanism of the Verwey transition in magnetite: Jahn-Teller distortion 

and charge ordering patterns. J. Phys.: Condens. Matter 18, 10427-10436 (2006). 
[17] Asamitsu, A., Tomioka, Y, Kuwahara, H. & Tokura, Y. Current switching of ressitive states in mag- 

netoresistive manganites. Nature 388, 50-52 (1997). 
[18] Sawa, A., Fujii, T., Kawasaki, M. & Tokura, Y. Hysteretic current-voltage characteristics and resis- 
tance switching at a rectifying Ti/Pro.rCao.sMnOa interface. Appl. Phys. Lett. 85, 4073-4075 (2004). 
[19] Yu, W. W., Falkner, J. C, Yavuz, C. T. & Colvin, V. L. Synthesis of monodisperse iron oxide nanocrys- 

tals by thermal decomposition of iron carboxylate salts. Chem. Commun. 2004, 2306-2307 (2004). 
[20] Zhou, Y, Jin, X. & Shvets, 1. V. Enhancement of the magnetization saturation in magnetite (100) 

epitaxial films by thermo-chemical treatment. J. Appl. Phys. 95, 7357-7359 (2004). 
[21] Perez-Diete, V. et al. Thermal decomposition of surfactant coatings on Co and Ni nanocrystals. Appl. 



Phys. Lett. 83, 5053-5055 (2003). 
[22] Zeng, H. etal. Magnetotransport of magnetite nanoparticle arrays. Phys. Rev. B 13, 020402(R) (2006). 
[23] Shvets, I. V. et al. Long-range charge order on the Fe3O4(001) surface. Phys. Rev. B 70, 155406 

(2004). 

[24] Burch, T. et al. Switching in magnetite: a thermally driven magnetic phase transition. Phys. Rev. Lett. 
23, 1444-1447 (1969). 

[25] Freud, P. J. & Hed, A. Z. Dynamics of the electric-field-induced conductivity transition in magnetite. 

Phys. Rev. Lett. 23, 1440-1443 (1969). 
[26] Duchene, J., Terraillon, M., Pailly, P. & Adam, G. Filamentary conduction in V02 coplanar thin-film 

devices. App/. Phys. Lett. 19, 115-117 (1971). 
[27] Gu, Q., Falk, A., Wu, J., Ouyang, L. & Park, H. Current-driven phase oscillation and domain-wall 

propagation in ^x^i-x^2 nanobeams. Nano Lett. 7, 363-366 (2007). 
[28] Salazar, A., Oleaga, A., Wiechec, A., Tamawski, Z. & Kozlowski, A. Thermal diffusivity of 

Fe3_a;Zna;04. IEEE Trans. Magn. 40, 2820-2822 (2007). 

II. ADDENDUM 

Supplementary Information is linked to the online version of the paper. 
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